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METHOD AND APPARATUS FOR 
ULTRASONIC CONTINUOUS, NON-INVASIVE 
BLOOD PRESSURE MONITORING 

BACKGROUND OF THE INVENTION 

This invention generally relates to methods and apparatus for 
continuous, non-invasive monitoring of blood pressure. 

A method for the determination of non-invasive beat-by-beat 
(continuous) systolic and diastolic blood pressure has long been desired in 
5 physiological monitoring. With this information, rapid changes in the 
physiological state of a patient can be much better managed. Automatic 
blood pressure cuffs can be used for this application, by inflating them as 
rapidly as is possible; however, this provides blood pressure data only at 1 to 
2 minute intervals, and each inflation can be painful for elderly or 
10 hypertensive patients. 

A very reliable technique for continuously measuring blood 
pressure is to insert a saline-filled catheter through the patient's vascular system 
to the point at which it is desired to perform the measurements. The catheter is 
connected to a pressure sensor, which measures the pressure in the vessel. An 
15 alternative method uses a catheter with a pressure sensor at the tip that directly 
senses the blood pressure. However, these techniques involve making an 
incision through the patient's skin and inserting the catheter into a blood vessel. 
As a consequence, this invasive procedure entails some risk of complications for 
the patient. 

20 An indirect, non-invasive process for continuously measuring blood 

pressure is based on the pulse transit time (PTT) which is the time required for a 
blood pressure pulse from the heart beat to propagate between two points in the 
vascular system. One apparatus for this technique includes an 
electrocardiograph that senses electrical signals in the heart to provide an 

25 indication when a blood pulse enters the aorta. A pulse oximeter is placed on an 
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index finger of the patient to detect when the blood pressure pulse reaches that 
location. The pulse transit time between the heart and the index finger is 
measured and calibrated to the existing blood pressure that is measured by 
another means, such as by the automated oscillometric method. Thereafter 
5 changes in the pulse transit time correspond to changes in the blood pressure. 
Generally, the faster the transit time the higher the blood pressure. Thus, 
changes in the pulse transit time can be equated to changes in the blood 
pressure. 

However, the electrocardiograph (ECG) senses electrical signals 
10 in the heart, which do not indicate the point in time when the blood pressure 
pulse actually leaves the heart upon the mechanical opening of the aortic valve. 
A time interval of varying length, known as the cardiac pre-ejection period (PEP), 
exists between peaks of the QRS wave of the electrocardiogram signal and the 
aortic valve opening. The inability of prior pulse transit time-based monitors to 
15 account for the cardiac pre-injection period resulted in an inaccurate 
measurement of the pulse transit time and thus blood pressure. 

In addition, changes in the compliance of the blood vessels also 
affect the pulse transit time. Chronic changes in arterial compliance occur due to 
aging, arteriosclerosis or hypertension. Arterial compliance can also change 
20 acutely due to neural, humoral, myogenic or other influences. Previous 
monitoring systems have been unable to separate changes due to compliance 
from changes due to blood pressure. As a consequence, some degree of 
inaccuracy has existed in calculating blood pressure from the variation of the 
pulse transit time. 

25 Some previous investigators in the area of continuous, non- 

invasive blood pressure (CNIBP) monitoring have favored tonometric 
methods. These methods do not require a blood pressure cuff, but do require 
some small mechanical device that applies pressure to an artery, along with 
some kind of vibration sensor for tonometric pressure estimation. Such a 
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device is described in published U.S. Patent Application Publ. No. 
2003/0149369 A1. Devices of this type have proven unreliable in practice. 

Still other previous investigators have used formulas derived 
from the Bramwell-Hill equation. These formulas relate blood pressure to 
5 measured arterial pulse wave velocity (PWV) and measured arterial blood 
volume. When the heart beats, it sends a pulse of pressure through the 
arterial system. This pulse propagates through the system by distending the 
elastic walls of the arteries, and this mechanism can be approximately 
represented by a linear wave equation. The basic expression for phase 
10 velocity of arterial pulse pressure waves is the same as that for the phase 
velocity of electrical waves propagating in a cable or transmission line, which 
has the same form of wave equation. The velocity is 



where v P is the phase velocity of the arterial pulse pressure wave. For cable 
15 waves, C is capacitance and L is inductance. For arterial pulse wave 
propagation, C is called the compliance of the artery, defined as the 
derivative of the lumen area with respect to the pressure, and is given by 

C = ^ = ^- (2) 
dP Eh v ' 

where A is the area of the arterial lumen, P is the pressure, r is the radius of 
20 the arterial lumen, h is the thickness of the arterial wall, and E is the modulus 
of elasticity of the arterial wall. This expression for C is derived in many 
standard textbooks, such as Hemodynamics, 2d edition, pages 96 - 98, by 
W.R. Milnor, published by Williams & Wilkins, 1989. L, which is the mass per 
unit length along the artery of the blood and represents the inertia of the 
25 blood that opposes the pressure pulse, is given by 
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L = -^ (3) 
7tr 



where p is the density of blood. If we substitute Eqs. (2) and (3) into Eq. (1), 
we get the famous Moens-Korteweg equation: 



If one substitutes Eq. (3) and the definition of C as a derivative into Eq. (1), 
we get what is called the Bramwell-Hill equation: 



p dA 



From this equation one can obtain a formula that linearly scales a measured 
arterial area into a pressure, with the slope based on measured PWV. If one 
10 puts all the pressure terms on one side and all the area terms on the other 
side (p is a constant) and then integrates, one gets: 



P(t)-P(0) = pv; In 



U(0). 



(6) 



Equation (6) is an approximation because the pulse wave velocity is a 
function of the arterial radius, and therefore the lumen area, and in the 
15 derivation of Eq. (6) it was assumed that v P is a constant with respect to A. It 
is also possible to use a linearized version of Eq. (6), given by: 

dv 2 

P(t)-P(0) = ^-(A(t)-A(0)) (7) 
A(0) 

In order to use Eq. (6) or (7) to generate P(t), the arterial pulse wave velocity 
(PWV) v P and the arterial (i.e., lumen) area A must be measured. In addition, 
20 the initial values P(0) and A(0) must be obtained, and the blood density, p, 
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must be replaced with a constant computed from cuff pressure data and 
ultrasound area and PWV measurements, all of which constitutes a 
calibration step. 

In previous methods, PWV measurements were typically taken 
by observing the pulse transit time between two widely separated sites, such 
as the heart and a fingertip. The pulse arrival times at the measurement sites 
are typically determined by impedance plethysmography or pulse oximetry. In 
one known prior method (disclosed in U.S. Patent No. 5,857,975), the time of 
the pressure pulse's origin at the heart is determined from an EKG signal. 
The required area measurements can also be obtained from 
plethysmography: volume measurements can be turned into areas by 
assuming a length. The initialization data can be obtained from a pressure 
cuff. The latter two measurements are not always used, however. 

There are also methods that use only pulse wave velocity; 
15 because these methods do not measure blood volume/area, they must use 
empirical relationships between PWV and blood pressure rather than Eq. (6) 
or (7). One such method is described in published U.S. patent application 
Ser. No. 2003/0167012. 

In the scheme described in U.S. Patent No. 5,857,975, 
20 measurement of area and calibration are replaced by use of an assortment of 
seemingly arbitrary factors. (This only highlights the fact that schemes that 
estimate blood pressure on the basis of PWV alone have no physical 
justification.) 

There are several grounds on which existing blood pressure 
25 calculations based on Eq. (6) might be criticized. One such criticism concerns 
the measurement of arterial blood volume. In particular, measurement of 
arterial area by plethysmography is confounded by the highly elastic nature 
of the veins. Impedance plethysmography measures total blood volume, and 
it is often applied to a limb such as the arm or leg. Since this measurement 
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includes the venous blood volume, it cannot be reliably processed to produce 
arterial lumen area. If the venous blood volume did not change, then it might 
be possible to calibrate for this effect, but the venous blood volume is 
strongly affected by the subjects position, since hydrostatic pressures can 
5 cause pooling of blood in the veins, which are more highly distensible than 
the arteries. 

Another criticism is that the PWV measurement resulting from 
measurement at widely separated sites is not directly applicable to Eq. (6). 
The PWV is being used in Eq. (6) to give information on the mechanical 

10 properties of the artery; however, the PWV depends not only on the elasticity 
of the arterial wall, but also on its thickness and the size of the lumen. Any 
PWV measurement that is made between widely separated sites is actually 
measuring PWV over a collection of branches of the arterial tree. Because 
these branches will have varying lumen area, the measured pulse transit time 

15 will reflect a composite of the component pulse wave velocities, and so it will 
not accurately reflect the mechanical properties of any particular arterial 
segment. It is possible to regard such a measurement as an approximation, 
but since Eq. (6) applies to a uniform, cylindrical tube, errors are inevitable, 
and their magnitude is hard to predict. 

20 Another criticism of a method based on Eq. (6) concerns the 

way in which the PWV measurements are made. At every bifurcation and 
change in lumen radius along the arterial tree, a pulse wave reflection will 
occur. These pulse wave reflections can change the apparent (measurable) 
pulse transit time. The PWV that is needed in Eq. (6) is that which would be 

25 observed in a very long, uniform tube, where there would be no reflections. 

A more fundamental criticism of blood pressure estimation 
based on Eq. (6) is that changes in physiological state can bring about 
changes in arterial wall elasticity. In general, this will change the measured 
PWV, and so the model of Eq. (6) will, to some extent, adapt to such 
30 changes. However, the operating point (P(0), A(0)) also depends on this 
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elasticity, as does any multiplicative constant used to replace p. If the elastic 
modulus of the wall changes, then the same pressure will be associated with 
a different area. If the operating point and calibration constant are obtained 
using a pressure cuff calibration, then any changes of arterial elasticity would 
5 in principle require re-calibration, and if no such re-calibration is done, then 
gross errors in pressure estimation can occur. For example, if the mean 
pressure goes down while the elasticity goes up (which means that the elastic 
modulus E goes down), then the mean area could also go up. In that case, 
without an adjustment of (P(0), A(0)), the mean of the estimated pressure 
10 would go up rather than down. 

There is a need for a method and means for acquiring data for 
use in the above-described blood pressure estimation scheme that 
overcomes the aforementioned drawbacks of prior art systems. 

BRIEF DESCRIPTION OF THE INVENTION 

The present invention overcomes the aforementioned drawbacks 
15 by using ultrasound to provide input data for the blood pressure estimation 
scheme. In a disclosed embodiment, a continuous, non-invasive blood pressure 
measurement is produced using an ultrasonic transducer array patch that is 
adhered or attached to the patient. The ultrasound transducer array is operated 
in an autonomous mode by a digital signal processor to obtain data from which 
20 blood pressure information can be derived. 

One aspect of the invention is a method for monitoring of blood 
pressure in an artery of a patient, comprising the following steps: (a) 
transmitting beams of ultrasonic wave energy that intersect the artery; (b) 
acquiring acoustic data by transducing ultrasound wave energy, transmitted 
25 and returned from the artery wall, the blood flowing through the artery and the 
tissue surrounding the artery, into electrical signals; (c) estimating the 
diameter of the artery based on a first set of the acquired acoustic data; (d) 
calculating the arterial lumen area based on the estimated arterial diameter; 
(e) estimating the velocity of a pulse wave traveling down the artery based on 
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a second set of the acquired acoustic data; and (f) computing the blood 
pressure as a function of the estimated lumen area, the estimated pulse wave 
velocity and respective reference values for blood pressure and lumen area. 

Another aspect of the invention is a system for autonomous 
5 monitoring of blood pressure in an artery of a patient, comprising an array of 
ultrasonic transducer elements, data processing means, and means for 
delivering signals derived from the output of the array to the data processing 
means, wherein the data processing means are programmed to perform the 
following steps: (a) controlling the array to transmit beams of ultrasonic wave 

10 energy; (b) beamforming acoustic data output from the array in response to 
impinging ultrasound wave energy transmitted and returned from the artery wall, 
the blood flowing through the artery and the tissue surrounding the artery; (c) 
estimating the diameter of the artery based on a first set of the acoustic data; (d) 
calculating the arterial lumen area based on the estimated arterial diameter; (e) 

15 estimating the velocity of a pulse wave traveling down the artery based on a 
second set of the acoustic data; and (f) computing the blood pressure as a 
function of the estimated lumen area, the estimated pulse wave velocity and 
respective reference values for blood pressure and lumen area. 

A further aspect of the invention is a method for estimating 
20 pulse wave velocity in an artery, comprising the following steps: (a) 
transmitting beams of ultrasonic wave energy that intersect the artery at two 
or more locations separated by a distance along the axis of the artery; (b) 
acquiring acoustic data by transducing ultrasound wave energy, transmitted 
and returned from the artery wall, the blood flowing through the artery and the 
25 tissue surrounding the artery, into electrical signals; and (c) estimating the 
velocity of a pulse wave traveling down the artery based at least partly on the 
acoustic data acquired at the first and second axial locations using an 
algorithm that corrects for pulse wave reflections. 
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Yet another aspect of the invention is a system for estimating 
pulse wave velocity in an artery, comprising an array of ultrasonic transducer 
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elements, data processing means, and means for delivering signals derived 
from the output of the array to the data processing means, wherein the data 
processing means are programmed to perform the following steps: (a) 
controlling the array to transmit beams of ultrasonic wave energy that 
5 intersect the artery at first and second locations separated by a distance 
along the axis of the artery; (b) beamforming acoustic data output from the 
array in response to impinging ultrasound wave energy transmitted and 
returned from the artery wall, the blood flowing through the artery and the 
tissue surrounding the artery; and (c) estimating the velocity of a pulse wave 
10 traveling down the artery based at least partly on the acoustic data acquired 
at the first and second locations using an algorithm that corrects for pulse 
wave reflections. 

Other aspects of the invention are disclosed and claimed below. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a block diagram of a CNIBP measurement system in 
1 5 accordance with one embodiment of the invention. 

FIG. 2 is a drawing showing a cross-sectional view of a typical 

cMUT cell. 

FIG. 3 is a drawing showing an isometric view of the cMUT cell 
shown in FIG. 2. 

20 FIG. 4 is a drawing showing an isometric view of a "daisy" 

subelement formed from seven hexagonal cMUT cells having their top and 
bottom electrodes respectively hard-wired together. 

FIG. 5 is a drawing showing a mosaic array comprising eight 
annular elements. 

25 FIG. 6 is a drawing showing the construction of a cMUT patch in 

accordance with one embodiment of the invention. In the depicted example, 
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the cMUT cells have been surface micromachined on the substrate. However, 
the cMUT cells could alternatively be bulk micromachined on the substrate, in 
which case the solid line at the interface of the cMUT and substrate layers 
could be changed to a dashed line. 

5 FIG. 7 is a drawing showing the CNIBP measurement concept in 

accordance with the disclosed embodiment of the invention. 

FIG. 8 is a drawing showing scanning of an ultrasound beam, 
which can be repeated frequently to generate a rectilinear two-dimensional 
image. 

10 FIG. 9 is a drawing showing steering of an ultrasound beam away 

from the normal. 

FIG. 10 is a drawing showing examples of matched filters for the 
anterior and posterior walls of the artery in accordance with one embodiment of 
the invention. 

15 FIG. 1 1 is a graph showing an example of the diameter estimation 

using the matched filter technique on a brachial artery. 

FIG. 12 is a drawing showing long-axis ultrasound data acquisition 
in accordance with an alternative method of diameter estimation. 

FIG. 13 is a graph of intensity versus range for a single frame of 
20 data acquired using long-axis ultrasound data acquisition. 

FIG. 14 is a graph showing an example of the diameter output 
using the peak tracking algorithm disclosed herein. 

Reference will now be made to the drawings in which similar 
elements in different drawings bear the same reference numerals. 
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DETAILED DESCRIPTION OF THE INVENTION 

In accordance with embodiments of the invention disclosed 
herein, ultrasound is used to provide input data for the blood pressure 
estimation scheme. The use of transcutaneous ultrasound can provide 
volumetric flow, arterial lumen area and pulse wave velocity information. In 
5 addition, ultrasound measurements can be taken in such a way that all the 
data describes a single, uniform arterial segment (e.g., the brachial and 
common carotid arteries can be used). This means that a computed area 
relates only to the arterial blood volume present. This also means that the 
measured PWV is directly related to the mechanical properties of the 

10 segment of elastic tube (artery) for which the blood volume is being 
measured. In the patient monitoring application, the operator of the 
ultrasound device is eliminated through the use of software that automatically 
locates the artery in the ultrasound data, e.g., using known edge detection 
techniques. Autonomous operation of the ultrasound system allows it to 

15 report blood pressure and blood flow traces to the clinical users without those 
users having to interpret an ultrasound image or operate an ultrasound 
imaging device. 

The use of ultrasound measurements in the continuous blood 
pressure monitoring application is enabled by micromachined ultrasonic 

20 transducer patch probe technology, which allows ultrasound data to be taken 
using a thin, lightweight probe that adheres to the patient's skin. Recently 
semiconductor processes have been used to manufacture ultrasonic 
transducers of a type known as micromachined ultrasonic transducers 
(MUTs), which may be of the capacitive (cMUT) or piezoelectric (pMUT) 

25 variety. cMUTs are tiny diaphragm-like devices with electrodes that convert 
the sound vibration of a received ultrasound signal into a modulated 
capacitance. For transmission the capacitive charge is modulated to vibrate 
the diaphragm of the device and thereby transmit a sound wave. 



11 



130897 



One advantage of MUTs is that they can be made using 
semiconductor fabrication processes, such as microfabrication processes 
grouped under the heading "micromachining , \ As explained in U.S. Patent No. 
6,359,367: 

5 Micromachining is the formation of microscopic structures using a 

combination or set of (A) Patterning tools (generally lithography 
such as projection-aligners or wafer-steppers), and (B) Deposition 
tools such as PVD (physical vapor deposition), CVD (chemical 
vapor deposition), LPCVD (low-pressure chemical vapor 

10 deposition), PECVD (plasma chemical vapor deposition), and (C) 

Etching tools such as wet-chemical etching, plasma-etching, ion- 
milling, sputter-etching or laser-etching. Micromachining is typically 
performed on substrates or wafers made of silicon, glass, sapphire 
or ceramic. Such substrates or wafers are generally very flat and 

15 smooth and have lateral dimensions in inches. They are usually 

processed as groups in cassettes as they travel from process tool 
to process tool. Each substrate can advantageously (but not 
necessarily) incorporate numerous copies of the product. There 
are two generic types of micromachining ... 1) Bulk 

20 micromachining wherein the wafer or substrate has large portions 

of its thickness sculptured, and 2) Surface micromachining wherein 
the sculpturing is generally limited to the surface, and particularly 
to thin deposited films on the surface. The micromachining 
definition used herein includes the use of conventional or known 

25 micromachinable materials including silicon, sapphire, glass 

materials of all types, polymers (such as polyimide), polysilicon, 
silicon nitride, silicon oxynitride, thin film metals such as aluminum 
alloys, copper alloys and tungsten, spin-on-glasses (SOGs), 
implantable or diffused dopants and grown films such as silicon 

30 oxides and nitrides. 

The same definition of micromachining is adopted herein. The systems 
resulting from such micromachining processes are typically referred to as 
"micromachined electro-mechanical systems" (MEMS). 

The use of a MUT patch allows the clinician to stick the 
35 transducer to the patient's skin. The MUT patch is lightweight and flat. For the 
purpose of illustration, a transducer patch will be described that is made up of 
capacitive micromachined ultrasonic transducers (cMUTs). However, it 
should be understood that the patch could instead employ pMUTs. An 
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embodiment will now be described that incorporates a MUT patch. However, 
it should be understood that the present invention encompasses not only a 
device, but also methods for continuous non-invasive blood pressure 
monitoring, and that the basic method is not limited in their application to 
5 devices that employ MUT technology. Instead, the basic method disclosed 
and broadly claimed herein can also be performed by appropriately 
programmed ultrasound imaging systems having probes comprising 
conventional piezoceramic transducer elements. 

FIG. 1 shows a simplified representation of a system for 
10 accomplishing the foregoing in accordance with one embodiment of the 
invention. A cMUT patch 30 is attached to the patient's skin with an adhesive 
and an inflatable cuff 34 is placed around one of the patient's arms. The 
patch is placed on the underside of the arm or in some other spot that gives it 
access to an artery that is relatively near the skin. The patch 30 is connected 
15 to a portable personal computer by a cable (not shown). The portable PC is 
small enough to be strapped to the patient. 

Based on acoustic data from the cMUT patch 30 and pressure 
data from the cuff 34, a digital signal processor (DSP) 32 [incorporated in the 
portable PC} autonomously computes various parameters, including the 

20 patient's estimated blood pressure. In addition, the DSP software forms and 
steers ultrasound beams over the volume of space in front of the cMUT patch 
30. Using both echo intensity and Doppler blood flow information, the 
software maps out the location of the vessel, its cross-sectional area, and the 
area increase at the time of passage of the systolic pressure wave. The pulse 

25 wave velocity can also be measured. This will yield useful information about 
the mechanical properties of the blood vessel wall. 

For this reason, the cuff 34 is periodically inflated by the DSP 
32 to provide the peak and steady state pressures at fixed points in time. 
More specifically, the DSP activates an air compressor 36 to inflate cuff 34. A 
30 pressure sensor 38 detects the pressure in the cuff 34. The analog output of 
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the pressure sensor 38 is converted to a digital value by an analog-to-digital 
converter 40, which digital pressure value is used by the DSP 32. Initially the 
DSP triggers cuff inflation in response to an activation command input by the 
system operator via an operator interface (not shown) of the portable PC. 
5 Thereafter, the DSP triggers re-inflation of the cuff whenever the estimated 
compliance value changes by more than a threshold. [The compliance is 
measured/estimated using arterial pulse wave velocity and mean arterial area 
using an algorithm that will be described later in detail.] The results of the 
computations are displayed on a display screen (not shown) of the portable 
10 PC. 

Although the pressure sensor and the transducer array in the 
disclosed embodiment are electrically coupled to the processor (or 
processors) by means of cables, obviously a wireless electrical coupling 
could be provided. 

15 Referring to FIG. 2, a typical cMUT transducer cell 2 is shown in 

cross section. An array of such cMUT transducer cells is typically fabricated on a 
substrate 4, such as a heavily doped silicon (hence, semiconductive) wafer. For 
each cMUT transducer cell, a thin membrane or diaphragm 8, which may be 
made of silicon nitride, is suspended above the substrate 4. The membrane 8 is 

20 supported on its periphery by an insulating support 6, which may be made of 
silicon oxide or silicon nitride. The cavity 16 between the membrane 8 and the 
substrate 4 may be air- or gas-filled or wholly or partially evacuated. A film or 
layer of conductive material, such as aluminum alloy or other suitable conductive 
material, forms an electrode 12 on the membrane 8, and another film or layer 

25 made of conductive material forms an electrode 10 on the substrate 4. 
Alternatively, the bottom electrode can be formed by appropriate doping of the 
semiconductive substrate 4. 

The two electrodes 10 and 12, separated by the cavity 16, form a 
capacitance. When an impinging acoustic signal causes the membrane 8 to 
30 vibrate, the variation in the capacitance can be detected using associated 
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electronics (not shown in FIG. 2), thereby transducing the acoustic signal into an 
electrical signal. Conversely, an AC signal applied to one of the electrodes will 
modulate the charge on the electrode, which in turn causes a modulation in the 
capacitive force between the electrodes, the latter causing the diaphragm to 
5 move and thereby transmit an acoustic signal. 

The individual cells can have round, rectangular, hexagonal, or 
other peripheral shapes. A cMUT cell having a hexagonal shape is shown in FIG. 
3. Hexagonal shapes provide dense packing of the cMUT cells of a transducer 
subelement. The cMUT cells can have different dimensions so that the 
10 transducer subelement will have composite characteristics of the different cell 
sizes, giving the transducer a broadband characteristic. 

Unfortunately, it is difficult to produce electronics that would 
allow individual control over such small cells. While in terms of the acoustical 
performance of the array as whole, the small cell size is excellent and leads 

15 to great flexibility, control is limited to larger structures. Grouping together 
multiple cells and connecting them electrically allows one to create a larger 
sub-element, which can have the individual control while maintaining the 
desired acoustical response. So a sub-element is a group of electrically 
connected cells that cannot be reconfigured. For the purpose of this 

20 disclosure, the sub-element is the smallest independently controlled 
acoustical unit. One can form rings or elements by connecting sub-elements 
together using a switching network. The elements can be reconfigured by 
changing the state of the switching network. However, sub-elements are 
hardwired into the design and cannot be reconfigured. 

25 For the purpose of illustration, FIG. 4 shows a "daisy" 

transducer sub-element 24 made up of seven hexagonal cMUT cells 2: a 
central cell surrounded by a ring of six cells, each cell in the ring being 
contiguous with a respective side of the central cell and the adjoining cells in 
the ring. The top electrodes 12 of each cell 2 are hard-wired together. In the 

30 case of a hexagonal array, six conductors 14 (shown in both FIGS. 3 and 4) 



130897 



radiate outward from the top electrode 12 and are respectively connected to 
the top electrodes of the neighboring cMUT cells (except in the case of cells 
on the periphery, which connect to three, not six, other cells). Similarly, the 
bottom electrodes 10 of each cell 2 are hard-wired together, forming a seven- 
5 times-larger capacitive transducer sub-element 24. 

Sub-elements of the type seen in FIG. 4 can be arranged to 
form a two-dimensional array on a semiconductive (e.g., silicon) substrate. 
These sub-elements can be reconfigured to form elements, such as annular 
rings, using a switching network. FIG. 5 is a drawing showing a mosaic array 
10 comprising eight annular elements. The drawing has been simplified by 
representing each subelement as a square, although it should be understood 
that the subelements could be of the type shown in FIG. 4. 

Reconfigurability using silicon-based ultrasound transducer 
sub-elements was described in U.S. Patent Application Ser. No. 10/383,990. 

15 One form of reconfigurability is the mosaic annular array, also described in 
that patent application. The mosaic annular array concept involves building 
annular elements by grouping sub-elements together using a reconfigurable 
electronic switching network. The goal is to reduce the number of 
beamforming channels, while maintaining image quality and improving slice 

20 thickness. To reduce system channels, the mosaic annular array makes use 
of the fact that for an unsteered beam, the delay contours on the surface of 
the underlying two-dimensional transducer array are circular. In other words, 
the iso-delay curves are annuli about the center of the beam. The circular 
symmetry of the delays leads to the obvious grouping of those sub-elements 

25 with common delays and thus the annular array is born. The reconfigurability 
can be used to step the beam along the larger underlying two-dimensional 
transducer array in order to form a scan or image. 

In accordance with one embodiment of the present invention 
shown in FIG. 6, an array of cMUT sub-elements is built on one silicon wafer and 
30 conventional metal oxide semiconductor (CMOS) switches and 
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preamplifier/buffer circuits are formed on a second silicon wafer to provide a 
cMUT patch having reconfigurable beamforming elements. An acoustic backing 
layer 18 is preferably sandwiched between the cMUT wafer 2/4 and the CMOS 
wafer 20 with vias 22 for passage of electrical connections between the wafers. 

5 The acoustic backing material 18 should have a composition 

that is acoustically matched to the cMUT substrate 4, to prevent reflection of 
the acoustic energy back into the device. In the case where the substrate is 
made of silicon, one example of a suitable backing material comprises a 
mixture of 96.3% (by mass) tungsten (of which 85% was 10 micron and 15% 

10 was 1 micron particle size) and 3.67% polyvinyl chloride (PVC) powders, as 
disclosed in U.S. Patent Application Ser. No. 10/248,022 entitled "Backing 
Material for Micromachined Ultrasonic Transducer Devices". The person 
skilled in the art will recognize that the composition of the acoustic backing 
material can be varied from the example given above, However, the acoustic 

15 impedance of the resulting backing material should be matched to that of 
substrate material. For example, if the substrate is silicon, the acoustic 
impedance should be approximately 19.8 MRayls ± 5%. 

The CMOS electronics preferably includes the transmit and 
receive circuits (including a respective transmit/receive switch for each cMUT 
20 sub-element) and at least a portion of the beamforming circuits. The CMOS 
electronics also include switches that enable reconfiguration of the sub-elements, 
allowing an aperture to be translated over the two-dimensional active area of the 
transducer. The shape of the apertures is determined by the desired steering 
angle for the ultrasound beam. 

25 For the autonomous blood pressure monitoring application, the 

artery will have to be identified without any intervention from an operator 
once the patch has been placed on the patient. The artery can be initially 
identified using the typical acquisition modes for color flow imaging available 
on current ultrasound machines. The color flow imaging is well suited for the 

30 task of highlighting the pulsating blood flow in a larger artery and can be 
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used along with the B-Mode images to provide an initial location of the artery 
for further tracking during the measurements. 

Once an artery has been identified, two types of data are acquired 
simultaneously using the cMUT patch, as shown in FIG. 7. The first data type is a 
5 sequence of B-mode images providing several short-axis views of the artery (in 
planes perpendicular to the direction of blood flow). Two or more short-axis views 
are taken along the length of the artery. The system processor can find the 
center of the artery automatically at each measurement location using edge 
detection techniques. The center of the artery can then be estimated by the 

10 system processor for any location along the length of the patch using 
interpolation. In addition to finding the artery center, the B-mode data is used to 
estimate the vessel diameter for the blood pressure calculation. More 
specifically, it is assumed that that the cross section of the lumen is a circle, so 
the estimated diameter can be used to calculate the lumen area A, which is then 

15 plugged into Eq. (6) or (7). 

The second data type is a set of M-mode lines to be used for 
Doppler processing. The M-mode lines are steered to intersect the center of the 
artery using information from the B-mode data. This beamsteering is done 
automatically by the system processor without operator intervention. For each M- 

20 mode line to be directed at locations along the artery that lie between B-mode 
planes, the beamsteering angle is adjusted by interpolating the B-mode arterial 
location data derived from those two B-mode image planes. The M-mode beams 
stay locked on the artery across the patch and in real time. The multiple M-mode 
lines can be used with no Doppler steering angle for determining the velocity of 

25 the vessel wall or with a Doppler steering angle for determining blood velocities. 
The M-mode data, either wall or blood velocities, will be used for estimating the 
pulse wave velocity needed for the blood pressure calculation. 

The two types of data are acquired by interleaving the acquisitions, 
possibly at different data rates, of the short-axis B-mode beams and the Doppler 
30 M-mode beams. An example of a repeating scan sequence for acquiring the data 
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could consist of a single B-mode beam at each location and all the M-mode 
beams. In the next pass through the sequence, the B-mode beams are translated 
across the patch while the M-mode beams remain in the same position. This is 
repeated until the desired width for the B-mode image is covered, generating a 
5 single frame of the B-mode data. Then the steps are repeated using any updated 
information about the vessel center. The result of the acquisition is M-mode data 
at a high repetition rate (~2 kHz), as required for Doppler processing, and the B- 
mode images at a slower rate (-100 Hz) determined by the required image size 
for automatic vessel tracking. 

10 The B-mode beams are translated across the cMUT patch by 

translating the annular array of activated sub-elements (seen in FIG. 5) across 
the patch as shown in FIG. 8. A uniform translation of the beamforming 
coefficients produces a new beam at a different location. Repeated frequently, 
this generates a rectilinear two-dimensional image. 

15 An additional bilinear term in the beamforming coefficients 

produces a beam directed away from the normal, as seen in FIG. 9. Such beams 
are transmitted in the M mode to acquire data for Doppler processing that 
computes the blood velocity based on the Doppler effect. 

M-mode beams that are transmitted generally normal to the artery 
20 are used to acquire data for Doppler processing of the arterial wall velocity 
during deflection of the wall as the pulse wave travels down the artery. The 
maximum deflection occurs when the wall velocity is zero. 

Blood pressure estimation is performed using a modified 
version of Eq. (7): 

25 PW-K, -^| - l]p. (0) - ^ V ^ (A. (t) - A(t)) (7') 
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In this formula, P(t) represents the estimated pressure at the time of the 
ultrasound measurements, v P (t) represents the estimated pulse wave velocity 
derived from the ultrasound measurements, and As(t) is the estimated average 
arterial area at systole from ultrasound measurements. A(t) can either be the 
instantaneous arterial area measured from ultrasound or the average area during 
some part of the heart cycle, such as diastole; in either case, P(t) is the 
corresponding blood pressure. The quantities identified with time zero are 
derived from the blood pressure cuff and from ultrasound measurements 
performed at the time of cuff inflation. The inflatable cuff measures systolic, 
diastolic and mean blood pressure over a period of 20 to 30 heartbeats. The 
measured systolic (maximum) blood pressure is P s (0) in Eq. (7'). Excursions of 
the lumen area A are detected using ultrasound to obtain the mean lumen area 
Am(0) in Eq. (7 1 ). Pulse wave velocity is estimated using methods described in 
detail below. The calibration blood pressure and the calibration lumen area and 
pulse wave velocity need not be measured concurrently. For example, the mean 
blood pressure can be measured over a span of 20 to 30 heartbeats and then 
the mean lumen area and pulse wave velocity can be measured over the next 20 
heartbeats. The value of K2 is computed so that the equation 



is satisfied for pressure and area measurements taken at nearly the same time. 

The term that multiplies P s (0) in Eq. (7') is there to modify the 
calibration time systolic pressure using the observed mean pulse wave velocity. 
This procedure is only valid during periods when the arterial tone does not 
change, and so this estimation scheme requires the use of the automatic re- 
calibration procedure described below. At the time of the first calibration, K1 is set 
to an empirically determined constant. When subsequent automatic re- 
calibrations occur, K1 is computed as 




(8) 
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K v p(°) , P,(T) 

1 v?(T)-v>(0) P s (0) 

where the time zero refers to the initial cuff calibration time and the time T refers 
to the subsequent cuff calibration time. 

The lumen area A(t) and the pulse wave velocity v p are measured 
5 to obtain the blood pressure P(t) using Eq. (7'). Measurements of the lumen area 
will typically be available at a far higher rate than measurements of the PWV. 

In order to estimate the instantaneous lumen area, it is assumed 
that the lumen is circular. The problem to solve is to determine the diameter of an 
artery in a short-axis ultrasound image sequence. In accordance with one 

1 0 embodiment of the invention, a two-dimensional matched filter method is used for 
estimating vessel diameter. In a short-axis ultrasound image, the lumen of the 
artery is black with a circular shape when the ultrasound scan plane is 
perpendicular to the direction of flow and becomes elliptical when the angle 
deviates from perpendicular. The edge that is seen in the ultrasound image is the 

15 boundary between the medial and adventitial layers of the artery. There is also a 
small portion of edge from the intimal layer, but this is only seen for a small arc 
near the anterior and posterior locations. 

To estimate the diameter from a circular image requires 
determination of two or more edge locations. If two locations are used, they must 

20 be constrained to be on a diameter. If three locations are used, the three 
parameters of a circle (center x and y location plus radius) can be found. If more 
than three locations are used, the circle parameter can be found to minimize an 
error. Similarly, to find the parameters of an ellipse requires four or more edge 
locations. To determine the diameter from an ellipse would require a correction 

25 for tilt of the scan plan from perpendicular. 

The embodiment described is a method for finding the edge 
locations in an ultrasound short-axis image to be used in determining the vessel 
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diameter. The approach is to construct a two-dimensional matched filter or 
template of the vessel wall at a certain location and then find the location in the 
image having the maximum correlation to the filter. The matched filter approach 
is an optimal method, in terms of signal-to-noise ratio, for finding a known signal 
5 in the presence of Gaussian white noise. The embodiment applies this concept 
to the artery wall detection. 

For the case of a circular-shaped vessel, the artery geometry is 
parameterized by an approximate diameter and edge thickness. A matched filter 
is designed to have the desired curvature for a portion of the edge and transition 

1 0 from a low value inside the lumen to a high value at the edge. The coefficients 
outside the desired edge are set to zero since this corresponds to the outer 
tissue, which varies significantly in intensity level. The pixel geometry in the 
image is also used in determining the coefficients of the two-dimensional 
matched filter. The approach will work for both linear probes that have 

1 5 rectangular pixels and for sector probes with wedge-shaped pixels. Examples of 
matched filters for the anterior and posterior walls of the artery are shown in FIG. 
10. 

Using the matched filter designed from the given vessel and pixel 
geometries, the filter is moved across the image and the correlation of the filter 

20 and the image is found at a set of pixel locations. The edge location is found by 
selecting the pixel with the highest correlation. Refinement of the edge location to 
sub-pixel accuracy (within a fraction of a pixel) can be performed by several 
methods. Interpolation could be done on the original image, thereby reducing the 
pixel size, or the pixel correlations could be fit with a surface and the location of 

25 the maximum on the correlation surface could be found more precisely. Another 
method would be to design a set of poly-phase filters corresponding to sub-pixel 
offset of the original filter. Correlation of these filters with the image at the current 
selected pixel could be used to determine the sub-pixel offset with the highest 
correlation. 
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The output of this method will determine a set of edge locations 
from which the diameter can be estimated. Since in most applications a 
sequence of ultrasound frames are acquired, the previous estimates of the 
diameter can be used to predict the diameter in the current frame and this value 
5 can be used to design or select the matched filter. In addition, the previous 
estimates of the edge locations can be used to predict the position of the edge in 
the current frame to define a small region of interest over which to search for the 
edge. An example of the diameter estimation using this technique on the brachial 
artery is shown in the FIG. 1 1 . 

Alternatively, a faster method for estimating the arterial area 
can be used. This alternative method uses only a single ultrasound line to 
measure the diameter and estimates the cross-sectional area from the single 
diameter measurement. This method is less reliable than the matched 
filtering method described above, but has the following advantages. (1) It is 
faster to collect the data for a single ultrasound line. (2) When performing the 
blood pressure estimation on an ordinary ultrasound device with a standard 
ultrasound probe, the short axis data is not available. 

This alternative method employs an algorithm that tracks arterial 
diameter in high-frame-rate long-axis ultrasound B-Mode data. Long-axis data is 
20 obtained by positioning the probe face in line with the artery, such that blood 
travels from one end of the probe face to the other. FIG. 12 is a diagram of a 
typical long-axis ultrasound data acquisition. 

In the case of CNIBP measurement, long-axis data acquisition is 
done at very high frame rates (-1000 frames per second, or fps). This high frame 
25 rate is needed for other necessary measurements, including pulse wave velocity. 
Since each beam fired requires time to travel the desired depth, reflect, and 
eventually die out, the number of beams fired and depth imaged will directly 
impact frame rate. For this reason, only four beams of B-mode ultrasound are 
fired to achieve -1000 fps. However, an image is not formed with this type of 
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acquisition; it is equivalent to multiple M-mode lines. FIG. 12 shows how these 
four beams are distributed: 

As shown in FIG. 12, an image is not formed with this high-frame- 
rate acquisition method, which eliminates the possibility of using adjacent beam 
5 information for noise suppression and filtering, as adjacent beams are separated 
by several millimeters. This is enough spatial distance for the artery to be in a 
different position and possibly have a slightly different diameter. 

The peak tracking algorithm described herein determines arterial 
diameter at every frame using a single beam of data. FIG. 13 shows a single 
10 frame of this data, B-Mode Intensity vs. Range (i.e., depth). The asterisks (*) in 
FIG. 13 denote the location of the anterior and posterior walls. At present, these 
wall locations are first selected manually as initial conditions. The algorithm 
tracks these locations in time, as described below, to determine arterial diameter 
for every frame. 

15 Peak Tracking Al gorithm 

As shown in FIG. 13, the data is inherently very noisy. The first 
step of the algorithm is to filter with a selectable-size Gaussian window, both 
spatially (along range) and temporally (along frames). The initial wall locations 
are identified manually (asterisk locations). Next, the algorithm fits a parabola to 

20 each of the initial wall locations and adjoining intensity values. The number of 
adjoining values used is selectable, but experimental results show a value of 3 
(for a total of 7 points; center ± 3) yields valid parabolas for most datasets. The 
range location of the maximum of each parabola is evaluated to obtain actual 
wall positions to sub-pixel accuracy. The difference of the anterior and posterior 

25 positions is the arterial diameter. 

For successive frames, the sub-pixel position is rounded to the 
nearest integer location. This is used as the next wall location estimate. It has 
been experimentally shown that the movement of these parabolic peaks allows 
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the algorithm to follow the arterial wall throughout the cardiac cycle. If the 
location of the parabolic peak is greater than 3 range samples away from the 
previous frame's estimate, or if the parabola is concave (indicating an invalid 
peak), the algorithm reports no movement and continues. This is done to prevent 
5 the algorithm from wandering from the actual wall location. 

The sub-pixel distance between anterior and posterior wall 
locations can be converted to absolute diameter in millimeters using knowledge 
of the range sampling for the acquisition. Time in seconds can be calculated in 
the same way using frame rate information. FIG. 14 shows a section of diameter 
1 0 output from the algorithm. 

More generally, the blood pressure estimation method should 
be able to adapt to physiological changes in arterial elasticity, which are part 
of the body's system for regulation of blood pressure. What is required is a 
method for sensing when a physiological change is the arterial elasticity has 
15 taken place and causing a re-calibration (that is, a cuff inflation). A PWV- 
based estimate of the arterial compliance C can be used as an indicator of 
physiological change. The Bramwell-Hill formula [Eq. (5)] gives a relationship 
between the PWV and the radius, parameterized by C, which is 

20 Since PWV and radius (or area) are measured, one can compute and track 
the estimated value of the compliance. The value of A in this formula is the 
mean area, and the value of K is computed using Eq. (9). The estimated 
compliance can be applied to the estimated variation in arterial area to 
estimate the pulse pressure, which is the difference between the systolic and 

25 diastolic pressures. When this estimated pulse pressure value changes by 
more than a threshold, a cuff inflation is triggered. 

25 



130897 



However, to compute the compliance used in conjunction with Eq. 
(9) to govern re-triggering of the blood pressure cuff, the actual pulse wave 
velocity, unbiased by reflections, is needed. As mentioned earlier, pulse wave 
reflections occur at every bifurcation and change in lumen radius. These pulse 
5 wave reflections can change the apparent (measurable) pulse transit time. In 
accordance with a further aspect of the present invention, a pulse wave velocity 
estimation technique is provided based on ultrasound measurements with 
corrections for reflections. 

In accordance with one embodiment of the invention, the DSP 

10 (or other processor) is programmed to perform PWV estimation with reflection 
correction. Processing steps for least-squares, frequency-domain estimation 
of arterial pulse wave velocity with model-based reflection correction will now 
be described. This method operates by measuring a disturbance caused by 
the passage of the arterial pulse pressure wave through an arterial segment. 

15 This disturbance could be the distension, or change in size, of the artery; it 
could be blood flow velocity; or it could be the velocity of the arterial wall. 
(The easiest technique is to measure wall velocity.) In any case, the 
measurement will be made at several sites along the length of the arterial 
segment of interest, and the measurements made at these sites will be 

20 processed jointly to obtain a pulse wave velocity estimate. Because the signal 
processing method is model-based, the disclosure will begin by describing a 
mathematical model that describes the observed disturbance in terms of the 
pulse pressure wave and its reflection. This description supplies sufficient 
background information to make the implementation of the scheme 

25 understandable. 

A Model nf th e Histurhanoe Components 

Because the sites of the ultrasound measurement lie along the 
length of the artery, the disturbance will appear to travel from one site to the 
next. Because the disturbance signal defined above responds to the pressure 
30 pulse, this signal propagates from site to site with the wave speed of the 
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pulse pressure wave. Thus, given the distance along the artery between 
measurement sites, the relative phases of the set of disturbance signals 
contain information about the pulse wave velocity. However, the direct 
measurement of this velocity is prevented by the presence of pulse wave 
5 reflections. If there were no pulse wave reflections, it would only be 
necessary to measure the delay of the disturbance as it propagates between 
measurement sites. With reflection, the shape of the disturbance changes 
between measurement sites, making delay estimation difficult. In order to 
estimate the pulse wave velocity in the presence of reflections, a model- 
10 based procedure is used that assumes that a linear reflection with unknown 
reflection coefficient occurs at a site an unknown distance distal from the 
most distal ultrasound measurement site. 

Suppose that the disturbance is observed at a site that is D 
millimeters upstream from a reflection site with reflection coefficient r. The 
15 total disturbance at that site can then be expressed as the sum of the forward 
and backward wave components: 

y( t ) = z(t) + rz(t-^) (10) 

where v P is defined as the pulse wave velocity. The term z(t) is the input to 
the arterial segment under consideration, and not, for example, the aortic 
20 pulse pressure waveform. Without knowledge of D and r t it would be difficult 
to resolve y(t) into forward and reverse components. 

Suppose that the disturbance is observed at N sites, separated 
by a distance d from one another, and that the most distal such site is D 
millimeters upstream from a reflection site with reflection coefficient T . The 
25 total disturbance at the i-th site can be expressed by an equation similar to 
Eq. (8). If all the components of the observed disturbances are expressed as 
delayed versions of the forward disturbance observed at the first site, then 
the disturbance at the n-th site is given by 
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for n = 1 to N. Equation (11) gives the n-th of N simultaneous equations, 
which all together may be expressed in matrix form: 
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5 This matrix expression allow one to solve for the three unknowns, D, r and 
v P . (The quantity d is known from the dimensions of the cMUT patch probe.) 
However, even though Eq. (12) represents N simultaneous equations in three 
scalar unknowns, there is still the unknown functional form of z(t) that stands 
in the way of its solution. This problem is resolved through a Fourier 

10 decomposition of z(t). Noting that t will be measured as a discrete time 
variable in the actual system, the discrete Fourier transform (DFT) is applied 
to the sampled version of z(t), transforming a finite segment of z(t) into a 
sequence of frequency components. (The DFT is computed by applying the 
well-known FFT algorithm.) After this has been done, a matrix equation 

15 similar to Eq. (12) is obtained for every coefficient in the computed DFT. 
Each of these matrix equations will represent N simultaneous equations in 
four unknowns, three real-valued and one complex-valued. Such an equation, 
whether underdetermined (for small N) or overdetermined (for large N), can 
be solved in the least squares sense to obtain an estimate of the pulse wave 

20 velocity. 

The observed signals y n (t), for n = 1 to N, are composed of 
linear combinations of delayed versions of z(t), which is a broadband 
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waveform. If one takes the DFT of a segment of a sampled version of y n (t), it 
will be possible to express it as a linear combination of the DFTs of z(t). 
Assuming that one has a segment of samples at times t = 0, T, .... mT, (M 
- 1)T, one can write the DFT as 



M-l 



-j 2 n km 
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(13) 



Substituting the definition of y n (t) from Eq. (10), one obtains 
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Now if one defines the DFT of z(mT) to be Z(k) f one has 
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Equation (14) is true for n = 1 N, and so can be written as a matrix 

equation, in the form of Eq. (12). Note that neither Z(k) nor the second 
exponential in the second term of Eq. (14) depends on n, so that the matrix 
equation is 



Y(k) = Z(k) s 



_2raik * 
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where e(O) is a vector defined as 
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1 



e j*(N-l) 

The parameter 0 is referred to herein as the phase increment of e, since it is 
the change in phase between successive components of e. The two phase 
increments given in Eq. (15) are negatives of one another, and are 
5 determined by the sample period, the distance between measurement sites, 
the DFT coefficient number, the number of samples used in the DFT and the 
pulse wave velocity. All of these quantities are known except the pulse wave 
velocity. 

The following six steps are executed by the system in order to 
10 estimate the pulse wave velocity. The first three steps produce the 
disturbance measurement, and the last three steps process the disturbance 
to produce a pulse wave velocity estimate. 

RtPp One- Data collPrtinn and Dopplfir shift estimation . 

First, an ultrasound data set is collected, consisting of a small 
15 number of l/Q lines, which pass through the artery of interest. These l/Q lines 
pass through the vessel at right angles to the arterial wall. For example, four 
or eight lines can be used. (Using more lines works better than using fewer 
lines.) This data set is taken repeatedly, at regular intervals, so that for each 
ultrasound line location, an ensemble of ultrasound lines suitable for color 
20 flow (mean Doppler shift) processing is available. The pulse repetition 
frequency for each such ultrasound line will be in the range of 1 to 5 kHz. The 
frequency of ultrasound used should be in the 8 to 12 MHz range. Since this 
data is taken repeatedly, an ensemble of any size may be obtained by simply 
storing in memory an appropriate number of the most recent ultrasound lines. 
25 This mode of operation is very similar to color M-mode, which is standard in 
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ultrasound imaging devices. An ensemble may be processed as each new 
ultrasound line becomes available, or at a lower rate. When Doppler shift 
estimates are based on ensembles that have significant data in common, the 
estimates are highly correlated. In the present invention, ensembles will be 
5 processed at the highest rate. 

Each of the ensembles of ultrasound l/Q lines generated by the 
data gathering process just described is processed to produce Doppler shift 
information related to the lateral (with respect to the blood flow direction) 
motion of the arterial walls. Each of the ensembles will first be filtered by a 

10 stationary tissue cancellation (high-pass) filter. This is a standard step in 
processing of ultrasound data to produce blood flow Doppler information; it is 
required for arterial wall motion processing because the increase in arterial 
diameter due to the pulse pressure wave can be quite small. Next, each 
ensemble will be processed to estimate its mean frequency using the 

15 standard autocorrelation method. This is also a standard step in color flow 
mapping. 

The result of the above processing is, for each ultrasound line, 
a set of Doppler shift estimates: one estimate for each range at which the 
ultrasound line is sampled. Range samples from the vicinity of the arterial 

20 wall closer to the ultrasound probe (the anterior wall) will feature an 
oscillating Doppler frequency indicating motion towards the probe. This is 
also true for much of the tissue between the probe and the artery. Range 
samples from the vicinity of the arterial wall farther from the ultrasound probe 
(the posterior wall) will produce Doppler data indicating an oscillatory motion 

25 away from the probe. 

Step Two- Segmentatio n of the Doppler shift data 

The next step in the processing is to segment the Doppler data 
into two non-overlapping sets of contiguous range samples, one indicating 
anterior wall motion and the other indicating posterior wall motion. Such a 



31 



130897 



segmentation can be accomplished by averaging the time samples 
associated with each range sample over some time window of fixed duration. 
Sites with an average motion toward the probe will be included in one 
segment and sites with an average motion away will form another. The 
5 duration of the time window over which the averaging takes place should be 
long enough to include at least one heart cycle: a period of 5 seconds 
duration would be long enough to encompass a single heart cycle at even the 
slowest heart rate. The segmentation will produce two sets of contiguous 
range samples, one set representing anterior wall motion and another 
10 representing posterior wall motion. 

Stpp Thrftfr Computation of thp observed disturbance 

s Once the segmentation has taken place, each new set of 
Doppler shift samples can be summed over each segment to produce two 
distinct quantities, one representing anterior wall motion and the other 
15 representing posterior wall motion. The difference between these two 
quantities 

z S (t)= IXk(0- 2> u (t) (17) 

k in posterior k in anterior 

segment segment 

will be called the measured disturbance for the i-th ultrasound line. In Eq. 
(17), Xi,k(t) represents the Doppler shift estimate computed for the i-th 

20 ultrasound line (out of four, say) at the k-th range (out of 128, say) at time t. 
In practice, t is quantized and represents the time of insonification. Zi(t) is an 
estimate of the time-varying velocity with which the anterior and posterior 
arterial walls are moving away from each other. This velocity is related to the 
rate of change of pressure within the vessel, and this changing pressure is 

25 caused by the passage of the arterial pulse pressure wave. 
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Step Four F requency decomposition of the disturbance observations 

Since z(t) is a random, almost-periodic signal that resembles a 
sequence of pulses, each period of the observed signal is Fourier 
transformed separately. Most of the information regarding the relative phase 
5 of the N signals will reside in the systolic response. During systole, the 
disturbance function will behave with the most energy. In order to segment 
the observations, the peak of one of them (the first ultrasound line, say) is 
detected and a set of samples is taken that includes some number of samples 
before and after this point. The resulting segment of recorded samples 

10 should represent a time segment whose duration is shorter than a heart 
cycle. For example, if the subject's pulse rate is 60 beats per minute and the 
sample rate of the ultrasound lines is 1000 insonifications per second, then it 
would be reasonable to choose 512 samples around peak systole to compute 
the DFT. These samples should be chosen so that there are a larger number 

15 of samples after peak systole than before, since the disturbance will be 
asymmetrical in that way. 

Note that all N of the disturbance functions must be Fourier 
transformed. The same set of time samples should be extracted from each for 
Fourier analysis. 

20 Step Five- Twn-step least-squares PWV solution for each DFT coefficient. 

If a vector of the N Fourier components with index k is formed, 
that vector will satisfy Eq. (15), according to the model. If the two unknown 
scalar coefficients of the two e(O) vectors in Eq. (15) are expressed as bi and 
b2, one can write Eq. (15) in the compact form 



25 Y(k) = b I e(* k (v p ))+ b 2 e(* k (-v p ))=E 



V 



(18) 



where the phase increment in the k-th DFT coefficient due to the pulse wave 
velocity has been written as <S> k(v p ) and the N-by-2 matrix E, whose columns 
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are the e(<t>) vectors, has been defined. Note that the matrix E is a function of 
the pulse wave velocity. 

Now Eq. (18) can be solved for [bi b2] T as follows: 



Lb 2 J 



= (e h e) 1 E H Y(k) (19) 



where the superscript H denotes Hermetian transposition. The squared error 
incurred by using the least-squares model to approximate the observation is 



Y(1c)-e(e h e)" 1 E h Y(k) 2 =|^I-E(E H E)T 1 E H jY(k) 



(20) 



which is a scalar function of the pulse wave velocity. The value of v P that 
minimizes Eq. (20) is the estimate of the pulse wave velocity that will be 
1 0 adopted. 

In order to compute the pulse wave velocity, values drawn from 
a range of possible PWVs from 1 to 10 meters per second are used to 
compute the observation error given by Eq. (20). This range of pulse wave 
velocities will be called the feasible range. (The Golden Section search is a 

15 well-known procedure that can be used to minimize the number of times the 
expression in Eq. (20) must be evaluated in order to find the minimizing 
value.) If the minimum observation error is associated with a value in the 
interior of the feasible range, that value is taken as the PWV estimate for the 
k-th DFT coefficient vector. If one of the two endpoints of the feasible range 

20 minimizes Eq. (20), then the estimate for that DFT index is ignored, since the 
true minimum is outside the feasible range. 

The operation described in this step is performed multiple times, 
once for each of a range of DFT coefficient indices. Only a relatively small 
number of indices, corresponding to low-frequency components, need be 
25 used for PWV computations. These indices correspond to frequencies at 
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which the power spectral density of the arterial pulse pressure wave is high. 
In practice, a 512-point DFT with coefficients indexed 3 through 18 has been 
used. 

Step Six: Combination of the PWV estimates for noisy observations . 

5 When the disturbance waveforms are simply delayed versions 

of the first-site observation with added reflections (this might be the case in a 
simulation), the technique above computes the true pulse wave velocity for 
each DFT index. In actual data, however, the measurements differ from each 
other due to added noise, the angle that the ultrasound beam makes with the 
10 artery and other factors. In this case, a different PWV estimate is obtained for 
each DFT index, in general. Since anywhere from eight to sixteen (or possibly 
more) such estimates will be available for each heart cycle, it will be 
necessary to have a method for combining the estimates into a single overall 
estimate for each heart cycle. 

15 In accordance with one embodiment, the PWV estimates are 

combined using what is called a trimmed mean. All the feasible estimates 
(those not at the endpoints of the feasible region) are sorted into numerical 
order. A number of the largest and the smallest are eliminated, and the 
remainder are averaged to get the PWV. The number eliminated depends on 

20 the number of feasible estimates observed. Alternatively, all the estimates 
from multiple heart cycles could be combined and processed in this manner 
to produce a single PWV estimate. 

The six steps above, executed repeatedly, produce a sequence 
of pulse wave velocity estimates at a maximum rate of one per heart cycle. 
25 These estimates may be used in the Bramwell-Hill equation, along with 
arterial area measurements, to estimate blood pressure and/or pulse 
pressure. In practice, it has been found that data must be combined over 20 
to 30 heart cycles in order to compute meaningful PWV estimates. 
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TimP rnrrpr.tions 

In the description given above, the equations were written as if 
all N ultrasound lines occurred at the same time. This is not the case. 
Actually, the insonifications occur in sequence, and there is a time skew 
5 between the data collection times. This skew can be accounted for by 
incorporating it as a fixed delay (for a given n) in Eq. (9). This value, if it is 
linear in the index n, can be accommodated by modifying the true value of d, 
the distance between the measurement sites, to one that accounts for the 
value of the ultrasound time skew. The remainder of the computations can 
10 proceed just as before, using this new value of the inter-measurement 
distance. 

While the invention has been described with reference to 
preferred embodiments, it will be understood by those skilled in the art that 
various changes may be made and equivalents may be substituted for 

15 elements thereof without departing from the scope of the invention. In 
addition, many modifications may be made to adapt a particular situation to 
the teachings of the invention without departing from the essential scope 
thereof. Therefore it is intended that the invention not be limited to the 
particular embodiment disclosed as the best mode contemplated for carrying 

20 out this invention, but that the invention will include all embodiments falling 
within the scope of the appended claims. 
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